Product Perfect Z^Z^-linear codes in Steganography 



Josep Rifa and Lorena Ronquillo 
Department of Information and Communications Engineering, 
Universitat Autonoma de Barcelona, 
08193-Cerdanyola del Valles, Spain. 
Email: Josep.Rifa@autonoma.edu, Lorena.Ronquillo@autonoma.edu 



Abstract — Product perfect codes have been proven to enhance 
the performance of the F5 steganographic method, whereas 
perfect Z2Z4-linear codes have been recently introduced as an 
efficient way to embed data, conforming to the ±1 -steganography. 
In this paper, we present two steganographic methods. On the 
one hand, a generalization of product perfect codes is made. On 
the other hand, this generalization is applied to perfect Z2Z4- 
linear codes. Finally, the performance of the proposed methods 
is evaluated and compared with those of the aforementioned 
schemes. 

I. Introduction 

Digital steganography is an information hiding application 
which consists of hiding data within a commonly used media 
in such a way that unintended recipients are not only unable 
to detect the presence of embedded data, but also given no 
reason for suspecting that anything is unusual. This is the main 
difference from encryption, which only prevents the adversary 
from decoding the message and not from suspecting that a 
secret message is being sent. 

Matrix encoding is a steganographic method introduced by 
Crandall [5] and analyzed by Bierbrauer et al. |2), which 
requires the sender and the recipient to agree in advance 
on a parity check matrix H, and the secret message is then 
extracted by the recipient as the syndrome (with respect to H) 
of the received cover object. This method was made popular 
by Westfeld fTO) , who incorporated a specific implementation 
using Hamming codes. The resulting method is known as the 
F5 algorithm, and it can embed t bits of message in 2* — 1 
cover symbols by changing, at most, one of them. 

The following two parameters are used to evaluate the 
performance of a steganographic method over a cover message 
of N symbols: the average distortion D = where R a is 
the expected number of changes over uniformly distributed 
messages; and the embedding rate E = which is the 
amount of bits that can be hidden in a cover message. Given 
two methods with the same embedding rate, the one with 
smaller average distortion is better. Following the terminology 
used by Fridrich et al. the tuple (D, E) will be called CI- 
rate. 

As Willems et al. in [ fTT| , we w dl also assume that a 
discrete source produces a sequence x = (x\, . . . , Xjsr), where 
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N is the block length, x t e N = {0,1,..., 2 B - 1}, and 
B E {8, 12, 16} depends on the kind of source (digital image, 
CD audio, etc). Let s G {1, . . . , M} be the message we want 
to hide into a host sequence x, which produces a composite 
sequence y = /(x, s), for y = (y 1} ... y N ) and y t € H. The 
sequence y is obtained from distorting x, and that distortion 
will be assumed to be of squared-error type (see fTT)). In these 
conditions, information can be carried by the least significant 
bit (LSB) or by the two least significant bits of each xi. An 
appropriate solution for the first case comes from applying the 
F5 algorithm [10], which has been improved in |7| by using 
the Kronecker product of the corresponding generator matrices 
of two binary perfect codes. The latter case is known as 
"±l-steganography" and the magnitude of changes is limited 
to 1, that is, yi = Xi + c, where c € {0, +1,— 1}. This 
case has usually involved the use of ternary codes (6j, (TTJ 
until the results from (8), which introduces a method based 
on perfect Z2Z4-linear codes. This kind of codes are not 
linear but have a representation using a parity check matrix 
that makes them as efficient as the Hamming codes. The 
steganographic method therein presented not only outperforms 
the one obtained by direct sum of ternary Hamming codes, but 
it also deals better with the extreme grayscale values problem. 
That is, the problem we may have when the steganographic 
method requires adding one unit to a grayscale value which 
already has the maximum allowed value 2 s — 1, or substracting 
one unit from a grayscale of value 0. 

Let H q (x) = Iog rm(ff(z) + x\og 2 (q - 1)) be the q- 
ary entropy function p| on the interval [0, (q — l)/q], where 
H(x) = — x\og 2 (x) — (l — x) log 2 (l — x) is the usual binary 
entropy function on the interval [0,1/2]. We call normalized 
embedding rate the ratio e = q D , where H q ~ (•) is the 
inverse of the g-ary entropy function H q (x). In the binary case, 
e will be computed by considering the binary entropy function 
H 2 (x), whereas in the ±l-steganography the ternary entropy 
function H 3 (x) is used. One of the purposes of steganographic 
methods is to approach the upper bound on the normalized 
embedding rate e subject to the constraint of an average 
distortion D. This upper bound on e for a fixed D is e < 1, 
and it is the same for any kind of steganography, be it binary 
or ±l-steganography. 

In this paper we propose a technique based on products of 
perfect Z^Z^-linear codes and compare its performance with 
that of the product binary perfect codes [7] and perfect Z2Z4- 



linear codes (8). 

The current paper has been organized as follows. Some 
basic concepts on perfect Z2Z4-linear codes, as well as the 
steganographic method based on these codes [8], are reviewed 
in Section [TTJ Then, Section [Til] reviews the product perfect 
codes method [7| and presents a generalization that enhances 
its performance. In Section IV this generalization is used for 
Z2Z4-additive perfect codes in ±l-steganography. Finally, the 
paper is concluded in Section [V] 

II. Perfect Z 2 Z 4 -linear codes and steganography 

Any non-empty subgroup C of ZS? x Z4 is a r Li7L^-additive 
code, where Z2 denotes the set of all binary vectors of length 
a and Z4 is the set of all quaternary vectors of length f3. Let 
be the usual Gray map from Z4 onto Z|, where 0(0) = 
(0,0), (f)(1) = (0,1), 0(2) = (1,1), and 0(3) = (1,0); and 
let <E> : Zj x Z4 — ^ Z2 be the extended Gray map given by 
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A Z2Z4-additive code C is isomorphic to an abelian struc- 
ture like 7L\ x Z\. Therefore, C has \C\ = 2~<A S codewords, 
where 2 1+s of them are of order two. We call such code C a 
Z2Z 4- additive code of type (a,/3;-f,8) and its binary image 
C = &(C) is a ^^-linear code of type (a, (3; 7, 6). Note that 
the Lee distance of a Z2Z4-additive code C coincides with the 
Hamming distance of the Z 2 Z 4 -linear code C, and that the 
binary code C may not be linear. 

The Z2Z4- additive dual code of C, denoted by C , is 
defined as the set of vectors in Zg x Z4 that are orthogonal 
to every codeword in C, being the definition of inner product 
in Zj x Z4 the following (see 
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where u, v G Z2 x Z4 and computations are made considering 
the zeros and ones in the a binary coordinates as quaternary 
zeros and ones, respectively. 

The binary code Cj_ = ^(C- 1 -), of length n = a + 2/3, is 
called the ^Z^dual code of C. 

A Z2Z4-additive code C is said to be perfect if code C ~ 
$(C) is a perfect Z 2 Z4-linear code, that is all vectors in Z^ are 
within distance one from a codeword and the distance between 
two codewords is, at least, 3. 

It is well known [4] that for any to > 2 and each 6 € 
{0, ■ ■ • , LtJ } tnere ex i sts a perfect Z2Z4-linear code C of 
binary length n = 2 m — 1, such that its Z2Z4-dual code is of 
type (a, (3; 7, 5), where a = 2™^ -1,0 = 2™- 1 - 2" 1 - 5 - 1 
and 7 = to — 2(5 (note that the binary length can be computed 
as n = a + 2f3). This allows us to write the parity check 
matrix He of any Z2Z4-additive perfect code C for a given 
value of 5. Matrix He can be represented by taking as columns 
all possible vectors in 1^ x Z 4 , up to sign changes. In this 
representation, there are a columns which correspond to the 
binary part of codewords in C, and (3 columns of order four 
which correspond to the quaternary part. We agree on a 



representation of the a binary coordinates as coordinates in 
{0, 2} G Z 4 . Let h,, for i G {1, . . . , a + f3}, denote the i-th 
column vector of He- 

Now we proceed to review how a perfect Z 2 Z4-linear code 
C = $(C) can be used in steganography. Consider its Z 2 Z4- 
dual, of type (a, (3; 7, 8), which gives us a parity check matrix 
He with 7 rows of order two and 5 rows of order four. 

Take iV = 2 m_1 and let x — (xi, . . . ,xn) be a source 
of grayscale symbols such that Xi G H = {0, 1, . . . , 2 B — 1}, 
where, for instance, B = 8 for grayscale images. 

We assume each grayscale symbol Xi is represented as a 
binary vector (f(B-i)i! ■ ■ • , v u, voi)> obtained by first repre- 
senting Xi in base 4 and then applying the Gray map to every 
quaternary symbol in the base 4 representation. For example, 
the grayscale value 239 is represented as the quaternary vector 
(3233), which then gives rise to the binary vector (10111010) 
after applying the Gray map 0. 

The TV-length packet x of symbols is translated into a 
vector w of a binary and j3 quaternary coordinates. The 
binary coordinates come from taking the least significant bit 
of the binary representation of Xy, that is Uoi> along with the 
two least significant bits Uij,i?oi of the following (a + l)/2 
grayscale symbols x^. The quaternary coordinates of w come 
from taking the two least significant bits of the last (3 symbols 
Xi and interpreting them as integer numbers _1 (fii,woi) m 
Z 4 . 

The obtained vector w G Zf x Zf is then distorted 
according to the matrix encoding method (5), | [T0[ in such 
a way that Hcv? T + e ■ = s holds, where s £ x Zj 
is the secret message we want to embed in x, the value of e 
can be {0,1,3}, the syndrome vector of w with respect to the 
parity check matrix He is Hcw T , and is a column vector 
in He- This method also deals with the extreme grayscale 
values in a rather efficient way: when a symbol Xi having an 
extreme value, be it or 2 B — 1, has to be distorted in a 
way that would lead its value out of the range defined by 
H, two other symbols are changed instead one magnitude. 
One of these symbols is always x-y. This method has C/-rate 

' 27V - 1 + — 1 + We refer the 
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reader to [8 1 for further details on this steganographic scheme. 

III. Product of perfect codes and steganography 

Let W q be a finite field of q elements, where q is a prime 
power. Let C be a Hamming code over F„ of length n = 

q m - 1 

and dimension n — to. Let Gc, He be, respectively, 
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a generator matrix and a parity check matrix for C. 

Definition 1: The Rronecker product of two matrices A = 
[a r ,t] and B = [bij] over F g is a new matrix A ® B obtained 
by changing any element a r t in A by the matrix a r t B. 

The g-ary code Cjj o ^ Ho , that is the dual of the code 
constructed by taking He ® He as generator matrix, is a 
[n 2 ,n 2 — to 2 ] code with covering radius p c ± = m 

- H C ®H C 

(see [9]). Codewords of Cjj ^ H can be seen as nxn matrices 



whose rows or columns are codewords in the q-ary Hamming 
code C. 

The q-ary code Cg g ®g c > constructed from the generator 
matrix Gc®Gc is a [n 2 , (n — m) 2 ] code with covering radius 
Pc Gc ®g c = n + 1 + 2(n — m — 1) J9|, whose codewords can 
be seen asnxn matrices where both rows and columns are 
codewords in the g-ary Hamming code C. We will refer to 
code Cg C 0G c as product code. 

There is an efficient steganographic method [7] which uses 
the above defined code Cg c ®g c > f° r 3 = 2, to embed data. 
Given two binary Hamming codes of the same length n = 
2 m — 1, their product is considered, which gives a linear code of 
length n = (2 m — l) 2 and dimension (n — m) 2 . Its codewords 
can be seen as (2 m — 1) x (2 m — 1) matrices where every 
row and every column are codewords in the binary Hamming 
code. The embedding scheme therefore consists of first taking 
blocks in the cover source of size (2 m — 1) X (2 m — 1), and 
then applying the F5 algorithm to every row and also to the 
first c columns, for 1 < c < 2 m_1 — 1. As proved in [7|, the 
performance one can obtain with this method is better than 
the one obtained by just using the conventional F5 algorithm 
on the corresponding codes with the same average distortion. 
We refer to [17] for further details on this method. 

Now, in this paper we will proceed with a generalization of 
the above procedure, by taking the product of more than two 
q-ary Hamming codes. 

Let C be a q-ary Hamming code of length n — q q Si > 
dimension n — m, with generator matrix Gc and parity check 
matrix He, as defined at the beginning of this section. Take 
the code C of all the nxn matrices such that all their rows, as 
well as their first column, are codewords in the q-ary Hamming 
code C. Code C is a [n 2 , n(n — m) — m] code, Cg c ®g c c 
C C Cft c ^ Hc with covering radius p c > = n + 1. 

For the sake of a well understanding, the following reason- 
ing will be limited to the binary case. However, a generaliza- 
tion to the q-ary case is straightforward. 

Just as the method based on the product of two Hamming 
codes from |7], this procedure consists of a row embedding 
and a column embedding steps. We will take the LSB bit of 
every grayscale symbol in the cover source and form blocks 
of size nxn, where n = 2 m — 1. Let Cij be the coordinate 
in the z-th row and j-th column of these blocks, where 
{l,...,n}. 

1) Rows Embedding: 

The matrix encoding standard procedure [5 1, [10] applied to 
every row lets us embed — bits with an average distortion 
of — coordinates, thus giving a Ci-rate of (— tt, ™)- 

n+l ' & o \ n+l ' n > 

2) Column Embedding: 

After processing all rows, we can embed ^ additional bits 
with an average distortion of ^rj by applying the same 
standard procedure to the first column. 
However, note that the following situations can happen 
when processing this column: 



• No coordinate needs to be changed in — =t- cases 

& n+l 

because the first column may already have, by chance, 
the desired value. 

• We may need to change a coordinate Cj,i in 
cases. In this case, the i-th row may have been already 
modified in the corresponding row embedding with 
a probability of ^xj, while it may have not been 
modified with a probability of ^j-p 

Let us consider the i-th row was modified in the 
coordinate qj, for j > 1. In this case, we will also 
have to restore the original value of Cjj and distort 
another appropiate coordinate c^t, for k € {2, . . . , n} 
and k ^ j, such that the distortion being introduced 
now by the column embedding is compensated and 
does not affect the embedding in the i-th row (see 
Lemma 2 from |7|). Note that this situation is also 
including the case in which the coordinate that was 
modified during the i-th row embedding is precisely 
the coordinate cj i we now need to change to embed 
data in the column. In summary, if the i-th row was 
modified, no matter in which coordinate, the column 
embedding step will introduce one distortion besides 
the ones introduced by the row embedding step. 

Otherwise, if during the column embedding we need 
to distort a coordinate a t \ and the i-th row was not 
modified, then we will also need to distort two more 
coordinates within the same row, Cjj and c^fc, for 
j,k G {2, ...,n} and j ^ k, to make up for this 
distortion. Hence, the column embedding step will be 
now introducing three changes. 

In short, we can leave invariant the average distortion of the 
row embedding step, but (g±§Mp+jj should be added for the 
embedding in the first column. Note that this is only a tight 
upper bound on the average distortion, as we will later show. 

By the method just described we can embed m bits into the 
first column and also in every row of the matrix; therefore, we 
embed (n + l)m bits in n 2 coordinates. The average distortion 
is upper bounded by ("+3)A^ l + 1 ) f or me coordinates in the first 
column and in each of the n rows. Summing this up, 

the average distortion is bounded by — — 5 

f , (»+3)/(n+l) \ 
n+l \ ~ n J' 

The method we propose in the present paper consists of 
repeating over and over the same procedure. Hence, we can 
generalize the computations of the average distortion and the 
embedding rate by using Gc 1 — Gc <8> {Gc ® • • ■ <8> Gc)- In 
each step Gc 1 = Gc ® Gc , on ly the first column in the 
first component Gc will be used to embed information. 

Let Di be the average distortion at the l-th step. As 
computed before, we have D\ = and D2 = ^xx(l + 

(n+3)/(»+l) y j n jjjg g enera J case we have; 

n + l 



n+3 



where f* — 

s n(n+l) 

Now, the overall average distortion can be computed as 
+ ■ + which converges asymptotically very 

fast to 
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As for the embedding rate, it can be computed as 



(l+n+n^H hn' 



which converges to 
• • • + n' _1 )m 



rt-l 



Finally, we obtain the asymptotical C/-rate 



m 
n — 1 

n 



n 2 — 3 ' n— 1 



Note that we are not able to generate an embedding scheme 
for any C/-rate but only for natural values of m. However, 
given any non-allowable parameter D for the average distor- 
tion, we can always take two codes with C/-rates (Di,Ei) 
and (£>2, E2), where D% < D < D2, such that their direct sum 
gives rise to a new CJ-rate (D, E), with D = AL>i + (l-A)D 2 
and E — \E\ + (1 — X)E 2 . 

A comparison of the normalized embedding rate e = 
q D y '- , where q = 2, as a function of the average distor- 
tion D for the introduced Kronecker product technique (KP- 
technique) and the standard matrix encoding procedures (5), 
| [T0) is shown in Fig.[T] As explained before, this plot has been 
made by first computing the allowable points (D, e), and then 
applying the direct sum between the codes corresponding to 
two contiguous points {D\, e\) and (D2, e.2), where D\ < D 2 . 

For the sake of simplicity, some particular cases which 
may produce a lower distortion have been omitted in the 
computation of the distortion in the above C/-rate. For this 
reason, the distortion D in that C7-rate is an upper bound. 
As an example of one of these cases, recall that the col- 
umn embedding step of our procedure is introducing three 
distortions when we need to change the c,.i coordinate and 
the i-th row was not modified in the row embedding step. 
Note, however, that there is one particular case in which we 
may need to introduce two distortions instead of three. This 
happens when there exist two other coordinates in the same 
column, c 3i i and Ck,i, for j,k £ {l,...,n} and j 7^ fc, 
whose distortion is equivalent to distort only c^\, and both 
j-th and fc-th rows were modified in the row embedding step. 
It is easy to see that we can distort coordinates Cj 1 and 
Ck,i instead of c^i, and perform afterwards the appropiate 
changes to compensate these distortions in the j-th and fc- 
th rows, respectively. Therefore, the column embedding step 
will be introducing two distortions besides the ones introduced 
in the row embedding step, and not three, as we previously 
stated. However, if no two other coordinates, cj t i and Ck,i can 
be found such that both j-th and fc-th rows were modified, 
then the column embedding step does actually introduce three 
distortions. We have implemented and executed a simulation 
of the embedding procedure described in this section which 
considers, among others, this particular case. For this reason 



the experimental results of the Kronecker product technique 
("KP-technique (simulation)" in Fig. [T} have lower average 
distortion than the results obtained from the above C/-rate 
(plotted as "KP-technique" in Fig. fTJ. 
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Fig. 1 . Normalized embedding rate e as a function of the average distortion 
D, of steganographic methods based on the matrix encoding procedure |5J, 
(TO), and on the Kronecker product technique, using an upper bound on the 
average distortion ("KP-technique") and using the experimental results ("KP- 
technique (simulation)"). 



IV. Product of perfect Z 2 Z 4 -linear codes 

The previous procedure deals with Hamming codes. Now, 
we will apply it to perfect Z 2 Z4-linear codes. Let C be a Z2Z4- 
additive perfect code of type (a, /3; 7, S) and binary length 
n = 2 rn — 1, for m > 2, and let Tic be its parity check matrix. 
Take the code C whose codewords are all the nx N matrices, 
where N = 2 m ~ 1 , such that all rows are codewords in C, and 
so is the first column after applying the inverse of the extended 
Gray map $. 

Take blocks of n x N grayscale symbols in the source, 



X n ,N 



where N = 2 m ~ 1 , and translate them into n vectors of a 
binary and /3 quaternary coordinates, as stated in Section [TT] 
and explained in depth in |8|. At the same time, the first 
coordinate of those n vectors is making up a binary vector 
of length n which can also be seen as a vector of a binary 
and f3 quaternary coordinates by means of the inverse of the 
extended Gray map $ (see Section [n}. Note that considering 
the n rows and the first column, we end up having n + 1 
different vectors of binary length n. 

The embedding procedure we will apply here is very similar 
to the KP-technique described in Section III Once the n + 1 



vectors have been translated into n + 1 vectors in Z2 x Zf , 
we will proceed by steps: first, we will apply the embedding 



scheme from (8j to every row, and then we will apply it 
to the first column. Each distortion in the rows will involve 
adding or subtracting one unit to/from a grayscale symbol, 
and this requires considering the possibility of having extreme 
grayscale values problems. Recall that, unlike vectors in rows, 
the vector in the first column is only made up of the least 
significant bit of n grayscale symbols and not of their two 
least significant bits. This means that any distortion over 
a coordinate in this vector will involve a flip in the least 
significant bit of a grayscale symbol Xi i, for i 6 {1, . . . , n}, 
which leads us to conclude that, unlike the rows embedding 
step, no extreme grayscale values problem will ever crop up 
during the column embedding step. 



of the average distortion and the embedding rate by taking 
the code whose codewords are all the /-dimensional matrices, 
where ! = nx(nx-"Xnx N), such that their rows and 
the first component of every dimension are codewords in the 
Z2Z4-additive perfect code C. 

Let Di be the average distortion at the l-th step. For 

2n+ VV 

the first steps we have D\ = jn+ip anc ' ^2 = 
' 1 ]. In the general case we have: 



Furthermore, as in the KP-technique from Section III dur- 
ing the column embedding step we have to consider different 
situations. The embedding method may require modifying a 
certain coordinate in the first column, and this coordinate 
may correspond to a row which was (or was not) modified 
during the row embedding step, or to a row that contains 
two distorted grayscale symbols, probably to deal with an 
extreme grayscale value problem. In any of these cases the 
action to be taken may vary, but still the aim is performing the 
appropriate changes in the affected row so that the distortion 
being introduced now by the column embedding step does 
not affect the embedding in the row. Take any two column 



(n+l) 2 



where £ 



n+3 



2n 
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(n + l) 2 
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n+3 



(n+^r)(«+l)' 
Now, the overall average distortion can be computed as 

, which converges asymptoti- 



cally very fast to 
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As for the embedding rate, it can be computed as 

( 1 +"+"'+;-+"'" 1)m , which converges to 1 ^ TT . 

Nn L 1 ' to N(n — 1) 

Finally, we obtain a C7-rate of [ 2 "*j^ ( ) , • mn 



N(n-1) 



h„h k of order four in matrix He, such that one Fi & sh ™s i a comparison of the normalized embedding 



vectors 

is the complementary of the other, that is hj = + 2, 
where 2 is the all-twos vector. The changes above mentioned 
will consist of considering that any distortion in coordinate 
Xi : i, for any i € {1, . . . , n}, can be compensated either by 
doing Xij-( a+ i)/ 2 + 1 and x iii _( Q , + i)/ 2 + 1 or by doing 
x ij-(a+i)/2 — 1 an d Xi,j-(a+i)/2 — 1- Note that, whenever 
possible, we will avoid modifying those grayscale symbols 
associated with column vectors in He that are complementary 
of themselves, because in this case we would have to distort 
the associated symbol in two units instead of one, which would 
not conform to ±l-steganography. 

By means of this method we can embed m bits into the 
first column x\ \, ■ ■ ■ , x n .\ and also in every row of the block. 
Since the first column is made up of n grayscale symbols 
and each row is made up of N symbols, we are actually 
embedding (n + l)m bits in nN = n(n + l)/2 symbols. 
It is easy to see that an upper bound of the average distortion 
for the symbols in the first column is ( -" + ^^ 1 " +1 - > . As for 
the symbols in each row, the average distortion is given by 

(see Section 



rate e 



H q - L {E) 
D 



, for q — 3, as a function of the average 
distortion D for the steganographic method based on perfect 
Z2Z4-linear codes (8), the one based on ternary Hamming 
codes (6), fTT) and the new method based on the product 
of perfect l-iL^-Xva&w codes. Recall that the distortion we 
have computed in the above C/-rate is an upper bound on 
the average distortion, meaning that lower distortion can be 
achieved in some particular cases, as it happened in the 
simulation results from Section [ill] 



2N 2 — (n+l) 2 

an upper bound for the average distortion is 
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In a similar way as we did in Section III we can repeat 
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Fig. 2. Normalized embedding rate e as a function of the average distortion 
D, of steganographic methods based on Z2Z4-additive perfect codes [8] 
("Z2Z4-add. perf. codes"), on ternary Hamming codes (6J, |1 1 1 and on the 
product of Z2Z4-additive perfect codes ("prod, of Z2Z4-add. perf. codes"). 



this method over and over and generalize the computations 



V. Conclusions 



The use of perfect TL-iL^-Xxxvtm codes in ±l-steganography 
was first proposed in (8). This method has a better perfor- 
mance compared to those based on the direct sum of ternary 
Hamming codes from (6) and (TTJ, and also deals with the 
extreme grayscale values more efficiently. 

In this paper we have presented a technique based on 
products of these perfect Z^Z^-linear codes. Therefore, the 
proposed method has all the advantages related to the per- 
formance and the processing of extreme grayscale values 
compared to the techniques based on the direct sum of ternary 
Hamming codes. Furthermore, we have shown that it performs 
better than the method in |8|. 
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